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Abstract 


An Euler shock-fitting marching code yields good agreement with 
semiempirically determined plume shapes , although the agreement de- 
creases somewhat with increasing nozzle angle and the attendant in- 
crease in the nonisentropic nature of the flow. Some calculations for a 
low-boom configuration with a sample engine indicated that , for flight at 
altitudes above 60000 ft, the plume effect is dominant. This negates the 
advantages of a low-boom design. At lower altitudes, plume effects are 
significant but of the order that can be incorporated into the low-boom 
design process. 


Introduction 

The plume that results from an underexpanded 
jet is a significant factor in assessing the feasibility 
of potential supersonic commercial aircraft because 
of its influence on the sonic boom signature and 
on aerodynamic performance. Consequently, reliable 
procedures for computing both the plume shape and 
the effects of the plume on the sonic boom signature 
are required. This report presents a comparison of 
plume shapes computed by several methods, together 
with an examination of the kinds of sonic boom 
effects the plume introduces. 

The primary reference in the current literature 
on plume shapes and plume sonic boom effects is 
reference 1. This reference describes a combined 
experimental-computational study to determine em- 
pirically the shapes of jet equivalent solid bodies, 
that is, an effective discrete jet free-stream slip line. 
These shapes were compared with plume shapes pre- 
dicted by two computational methods. Reference 1 
also includes some sample calculations to illustrate 
the plume sonic boom effect. 

The present report describes calculations that 
yield a closer correlation with the semiempirically de- 
termined shapes than the methods described in ref- 
erence 1. Also included is a study of the influence 
of flight altitude and configuration geometry on the 
problem of plume sonic boom effects with some illus- 
trative examples. 

Symbols 

A{x) equivalent area distribution 

D external nozzle exit diameter, d + Lip 

thickness 

d internal nozzle diameter at jet exit, 2 h 

F value of Whitham F-function 

h internal nozzle radius at jet exit 


L airplane reference length 

M Mach number 

p static pressure 

Ap = p- Po 

r radial coordinate 

t time 

t = 0 time at which a pulse, traveling at 
ambient sound speed, arrives at 
ground 

x coordinate in direction of nozzle axis 

a nozzle expansion angle 

j3 nozzle boattail angle 

e dummy integration variable 

Subscripts: 

j jet 

o undisturbed condition 

oo free stream 

Plume Calculations 

The plume shapes were computed with the Euler 
shock-fitting marching code described in reference 2. 
Calculations were made for some, but not all, of 
the nozzle shapes and conditions included in the 
investigation described in reference 1. Calculations 
were not made for the boattailed shapes because 
the plume code could not treat this case without 
modification. 

The type of nozzle that was tested is depicted in 
figure 1. The nozzle divergence angles a for the non- 
boattailed shapes were 7.28°, 9.06°, and 11.50°. The 
measurements of reference 1 indicated that, for the 
two larger divergence angles, an internal shock was 
generated inside the nozzle. This shock is apparently 
associated with the curvature discontinuity that oc- 
curs where the circular-arc throat section meets the 
conical nozzle section. The experimental results in- 
dicated that this shock was relatively weak and its 


presence was not considered in the direct calculation 
of the plume shape or in the calculation of the exit 
Mach number. 

Figure 2 shows a set of calculations for plumes 
emitted from the nozzle with a — 11.5°. In this fig- 
ure, the solid line indicates the semiempirically de- 
termined jet equivalent plume shape. It was obtained 
by measuring the flow conditions along a line outside 
the plume and parallel to the nozzle axis. Then, with 
these initial conditions, a characteristic net was com- 
puted inward, taking the streamline that matches the 
nozzle lip as the jet equivalent solid body. 


No further attempt was made to investigate the 
various sources of error. The accuracy of the com- 
puted shapes is probably consistent with the accu- 
racy of the approximate procedures that arc used 
to obtain the effective area distributions that are 
required for sonic boom calculations (ref. 4). For 
the purpose of comparing computed shapes with the 
semiempirical shapes of reference 1, a simple adjust- 
ment formula could be applied, since the discrepancy 
increases with nozzle divergence angle in a systematic 
manner. 

Sonic Boom Considerations 


The dash lines in figure 2 denote calculations, 
reported in reference 1, by a linearized technique 
described in reference 3. The long-dash -short-dash 
lines, also from reference 1, denote calculations by 
a method of characteristics assuming straight shock 
lines and conical nozzle flow. The circles denote 
calculations by the Euler method of reference 2. 

All these calculations predict plume shapes larger 
than the semiempirically determined shapes. The 
discrepancy is significantly smaller for the Euler 
calculation. 

Figure 3 shows a similar comparison of computed 
and semiempirically determined shapes for a nozzle 
w r ith a = 9.0G°. For this nozzle, the points calculated 
by the Euler method represent a close approximation 
to the empirical shape, although it is still slightly 
overpredicted. 

The nozzle with a = 7.28° was tested at only one 
pressure ratio. For this case, shown in figure 4, the 
Euler method yields an excellent representation of 
the empirically determined shape. 

Thus these results indicate that as the nozzle ex- 
pansion angle increases, the correlation of the Euler 
calculations and the empirical results tends to de- 
teriorate. As the nozzle angle increases, the obser- 
vation was made that an internal shock forms. If 
this shock significantly reduces the Mach number at 
the nozzle lip, this effect would account for at least 
part of the disagreement, since the exit Mach num- 
ber is computed in reference 1 with the assumption 
of isentropic nozzle flow. However, a number of other 
factors could be involved. These factors include vis- 
cous effects, the effect of the bluntness of the nozzle 
lip, and assumptions involved in initializing the Euler 
plume code. Some attempts were made to determine 
the precise nozzle exit flow by computing the inter- 
nal nozzle flow field by time-relaxation Euler codes, 
but these attempts failed to yield the type of shock 
structure observed in the experimental results. 


Once the jet plumes have been modelled as equiv- 
alent solid bodies their sonic boom effects can readily 
be determined- The equivalent solid body is simply 
treated as an extension of the engine nacelle. Then 
an equivalent area distribution obtained from pro- 
jected Mach slice cuts through the configuration is 
generated. This area distribution is added to an 
equivalent area distribution due to lift to obtain a 
total area distribution A(x). This area distribution 
is related to the sonic boom wave shape through the 
“F- function” which is defined by the formula 


F(x) = 


f x A"{x - e) dc 
27T J 0 yjx — t 


The wave shape is determined directly from the 
F-function, and the ground level signature is de- 
termined by a propagation and shock-fitting code. 
These concepts arc explained in greater detail in the 
standard sonic boom literature. (Sec, for example, 
ref. 5.) 

In view of the wide variety of approaches to 
configuration geometry and aircraft cruise conditions 
that are being considered for civil supersonic aircraft, 
it would not be practical at this point to attempt to 
compute plume effects for each case. However, a few 
calculations in addition to those already reported in 
reference 1 may serve to indicate the types and orders 
of magnitude of the effects to be expected. 

Some factors involved are the flight Mach num- 
ber and altitude, the configuration geometry, and the 
manner in which the engine nacelles are integrated 
into the geometry. Configurations that are designed 
for diminished sonic boom levels are usually laid out 
in such a way that the lift is distributed longitudi- 
nally to the greatest possible extent. The result is 
a planform somewhat like that shown schematically 
in figure 5(a). This geometry may be compared with 
that of a conventional supersonic fighter-type design 
with a more concentrated lift distribution (fig. 5(b)) 
or with a conventional civil transport configuration 
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without highly distributed lift but with an exten- 
sively notched arrow wing (fig. 5(c)). 

For a relatively high flight altitude, the ambient 
pressure is so low that the nozzle flow is greatly un- 
derexpanded and consequently large plume effects 
are realized. For example, figure 6 shows the plume 
shape computed for a high performance afterburn- 
ing engine with assumed ambient pressure equiva- 
lent to that at 60 000 ft (A/qc — 3.0, pj/poo = 8.86). 
Flow parameters for this engine were computed by 
the method of reference 6. The pluming is signifi- 
cant, and consequently, the plumes for the four en- 
gines would contribute a significant equivalent area. 
The actual effect on the F-function and ground level 
signature is shown in figure 7 for a low-boom con- 
figuration of the type shown in figure 5(a). Fig- 
ure 7(a) represents the F- function and signature with 
a cylindrical afterbody (no pluming) assumed, and 
figure 7(b) shows the corresponding results with the 
plume of figure 6. The plume has a dominant effect 
on the sonic boom signature for this case. The com- 
pression associated with the plume causes a shock of 
such magnitude that it moves forward and overrides 
the nose shock. The ground level overpressure is in- 
creased from a level slightly under 1 psf (fig. 7(a)) to 
about 2.2 psf (fig. 7(b)). It would not be feasible to 
attempt to tailor the configuration to allow for the 
additional area associated with this plume, since this 
area is so large. 

However, the problem is mitigated at lower alti- 
tudes where the ambient pressure is higher and the 
plume is consequently smaller. Figure 8 shows the 
plume shape (with the same engine) for flight at 
55 000 ft and M 0 0 = 2.1. This plume is considerably 
smaller than that shown in figure 6, but it still has 
a significant effect on the F-function and signature. 
Figure 9(a) gives these results for the configuration 
with no plume, and figure 9(b) gives the correspond- 
ing results with the plume effect included. In this 
case the plume compression significantly alters the 
F-function and the signature, but the effect is much 
smaller than in the previous case. 

A cursory effort was made to reduce the plume 
effect further by staggering the engine nacelles a dis- 
tance of 10 ft and making modest changes (about 
20-percent local variation) in the fuselage area dis- 
tribution. The result, shown in figure 10, is some 
reduction of the plume effect. The maximum over- 
pressure is reduced from 1.8 to 1.6 psf. Further re- 
duction could probably be realized through a more 
systematic design approach (ref. 7), which might in- 
corporate techniques such as those of references 8 
and 9. 


Reducing the flight altitude from 65 000 to 
55000 ft is so effective in reducing the plume effect 
that it would appear that, by extrapolation, a further 
reduction in cruise altitude, say to 45 000 ft, would 
render the plume effect negligible. However, such 
an extrapolation cannot be made. The high perfor- 
mance afterburning engines are not appropriate for 
the lower altitude flight, which would be associated 
with a lower Mach number. 

Consequently, a sample case was computed for 
M =1.6 flight at = 1.6 and 45 000 ft with con- 
ventional turbo engines having an internal nozzle ex- 
pansion of about 58 percent. The results are shown 
in figures 11 and 12. Figure 11 shows the plume 
shape. Figure 12(a) shows the F-function and signa- 
ture with the assumption of a cylindrical plume, and 
figure 12(b) gives the corresponding results with the 
plume shape of figure 11. Again, a slight reduction 
in the plume effect can be realized by staggering the 
engines and tailoring the fuselage (fig. 13). This tai- 
loring is substantial, about a 40-percent variation in 
cross-sectional area, as is shown in figure 14. 

The results given in figures 9 and 1 1 may be com- 
pared with some calculations shown in figure 20 of 
reference 1. Those calculations demonstrate that, for 
configurations like the ones depicted in figures 5(b) 
and 5(c), relatively small jet plumes actually have a 
favorable effect on the sonic boom signature. The 
explanation for this difference involves both the am- 
plitude of the aft expansion region of the F-function 
and the location of the plume effect relative to this 
expansion region. Low-boom configurations have a 
relatively gradual, low amplitude expansion (as illus- 
trated in figs. 7(a), 9(a), and 12(a)), which results 
in a ground level amplitude of about 1.0 to 1.5 psf. 
When the compressive effect of the plume is imposed 
on the expanding flow, it is only partially cancelled 
by the expansion and leaves a secondary compres- 
sion. On the other hand, a conventional configura- 
tion design, like that of figure 5(b) or (c), has a more 
sudden and larger amplitude expansion region in the 
F-function; this leads to a ground level amplitude on 
the order of 2.5 psf (ref. 1). In this case, the plume 
compression can be completely submersed in this 
expansion. 

Figure 5(b) also illustrates another means that 
can be effective in controlling the location of the 
plume. Engine nacelles mounted on the fuselage can 
be situated in the optimum longitudinal position. 

It should be emphasized that the configurations of 
figures 5(b) and (c) are not low-boom designs. Even 
with some reduction of the tail wave strength, the 
concentrated lift distributions associated with such 
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designs generally yield unacceptably high compres- 
sive overpressures. 

Concluding Remarks 

An Euler shock-fitting marching code yielded rel- 
atively good agreement with semiempirically deter- 
mined plume shapes, although the agreement de- 
creased somewhat with increasing nozzle expansion 
angle. Some evidence indicates that the discrepancy 
may be attributable to nonisentropic internal noz- 
zle flow which is not accounted for in initializing the 
plume calculation. 

Some calculations were carried out to obtain a 
general assessment of the nature and orders of mag- 
nitudes of the plume effects on the sonic boom 
signature. The calculations were for a low-boom 
configuration with a high performance engine. The 
results indicated that, for flight at altitudes above 
60 000 ft, the plume effects were dominant; but for al- 
titudes below 55 000 ft, they were significant but not 
dominant. 

Some factors associated with incorporating the 
plume shape into the configuration design were also 
discussed. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
January 24, 1992 
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Figure 1. Nozzle geometry. 
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Semiempirical solid body Method of characteristics 

Inviscid jet boundary (ref. 3) O Euler shock fitting 



Semiempirical solid body 
Inviscid jet boundary (ref. 3) 
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Figure 8, Computed plume 




(a) Cylindrical jet. 

Figure 9. F-functions and ground level signatures for low-boom configuration with cylindrical jet and computed 
plume of figure 6 at M 0 0 = 2.1 and altitude of 55000 ft. 
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11. Computed plume shape for 



t 

(a) Cylindrical jet. 


Figure 12. F-functions and ground level signatures for low-boom configuration with cylindrical plume and 
computed plume of figure 11 at Moo = 1-6 and altitude of 45000 ft with turbo engine. 



Figure 13. F-function and signature for configuration with computed plume of figure 11 at M a c 
altitude of 45 000 ft with fuselage modified and engines staggered. 
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